Impact of strange quark matter nuggets on pycnonuclear reaction rates in the crusts 
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This paper presents an investigation into the pycnonuclear reaction rates in dense crustal matter 
of neutron stars contaminated with strange quark matter nuggets. The presence of such nuggets 
in the crustal matter of neutron stars would be a natural consequence if Witten's strange quark 
matter hypothesis is correct. The methodology presented in this paper is a recreation of a recent 
representation of nuclear force interactions embedded within pycnonuclear reaction processes. The 
study then extends the methodology to incorporate distinctive theoretical characteristics of strange 
quark matter nuggets, like their low charge-per-baryon ratio, and then assesses their effects on the 
pycnonuclear reaction rates. Particular emphasis is put on the impact of color superconductivity 
on the reaction rates. Depending on whether or not quark nuggets are in this novel state of matter, 
their electric charge properties vary drastically which turns out to have a dramatic effect on the 
pycnonuclear reaction rates. Future nuclear fusion network calculations may thus have the potential 
to shed light on the existence of strange quark matter nuggets and on whether or not they are in a 
color superconducting state, as suggested by QCD. 
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I. INTRODUCTION 



First proposed in 1934 and subsequently observed in 
1967, neutron stars are the remnants of supernova ex- 
plosions of aging supergiant stars. They are extremely 
massive at approximately 1.4 times the mass of the sun. 
The interior structure of a neutron star is still the sub- 
ject of debate, as neutron stars are tremendously compact 
objects, crushing their mass into a radius of roughly ten 
kilometers. This makes them some of the densest objects 
known to man, with core densities possibly more than ten 
times the density of atomic nuclei. This extreme com- 
pression provides a high-pressure environment in which 
numerous subatomic particle processes are believed to 
compete with each other and novel phases of matter may 
exist. The most spectacular ones stretch from the gener- 
ation of hyperons and baryon resonances P,[3, to quark 



deconfinement [1, 
boson condensates TT,!! 

[l, 2, 16, 17, 18, 19, 2Q1.) It has also been suggested that 
strange quark matter |2l|, [22| , made up of roughly equal 



10 1 , to the formation of 
15|. (For overviews, see 



uarks, may be more 



26|,|27|,|28|. This in- 



numbers of up, down and strange qi 
stable than ordinary atomic nuclei 
triguing possibility is known as the strange quark matter 
hypothesis. In the latter event, neutron stars could be 
made almost entirely of strange quark matter rather than 
confined hadronic matter [2^,13, [1^ . If quark matter ex- 
ists in neutron stars it ought to be a color superconductor 
[2^ [sol , [sil . [3^ . This fascinating possibility has renewed 
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tremendous interest in the physics and astrophysics of 
quark matter [13, [13, [H, [Ml, [ll . 

The idea that quark matter may exist in neutron stars 
is not new [1 i, E H, [3, H, ll]. With the dissolution of 
nuclei in the cores of neutron stars, the crush of gravity 
may allow even protons and neutrons to be broken into 
their constituent components. The baryons at this point 
exist in such close proximity that their quarks are effec- 
tively free. With up and down quarks being the least 
massive, they would be the first to appear in this de- 
confined state. Furthermore, given the relatively small 
mass of the strange quark, it is possible that highly ener- 
getic up and down quarks from nucleons may convert 
to low-energy strange quarks at about the same den- 
sity This mixture of up, down, and strange quarks 
is referred to as strange quark matter if its energy per 
baryon is less than that of ordinary nuclear matter, as 
expressed by the strange quark matter hypothesis. An- 
other distinguishing feature between nuclear matter and 
strange matter is that the latter is self-bound and elec- 
trically charge neutral. For that reason strange mat- 
ter objects can exist stably for a tremendous range of 
baryon numbers, A, ranging from ^ 10^ to 10^^. The 
low baryon-number end is formed by strangelets (strange 
nuggets) dH, m, m, m, 1131 , whUe the high baryon- 
number end is populated by strange (quark matter) stars 
[Ei,[I3,|2l[2l[2l. 

If strange quark matter exists, it may be present in any 
of the neutron star regions - in the outer and inner crust 
through accretion from a companion compact stellar ob- 
ject and in the core through state changes. The base of 
the outer crust is characterized by a transition known 
as neutron drip, where the attractive nuclear force that 
holds subatomic particles together in nuclei is saturated 
and neutrons subsequently leak out of the nuclei. The in- 
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ner crust is a composite of these free neutrons and a lat- 
tice of neutron-rich matter, making this region a suitable 
environment for the study of high-density nuclear interac- 
tions, known as pycnonuclear reactions when a lattice is 
involved. The lattice structure is a critical feature, as the 
extreme gravitational force of neutron stars can offset the 
kinetic and potential energies of the constituent matter 
in the crust of neutron stars, with the ground state lattice 
configuration balancing the competing forces. Within 
this framework, pycnonuclear reactions are distinctive in 
that they can occur at very low temperatures. In the in- 
ner crust, it is possible that pycnonuclear reaction rates 
between SQM and normal atomic matter may be very 
different from those between ordinary nuclei. Along this 
line of inquiry, this is a first step in the examination of py- 
cnonuclear reaction rates between normal nuclear matter 
and a simple model incorporating some distinct charac- 
teristics of strange quark matter. 

To this end, the current undertaking recreates re- 
cent efforts to quantify pycnonuclear reaction rates, ef- 
forts that included ion-ion interaction models system- 
atically constructed from nucleon-nucleon interactions. 
This study then widens the scope, allowing for the inclu- 
sion of theoretical characteristics of strange quark mat- 
ter. The final encoded methodology calculates rudimen- 
tary pycnonuclear reaction rates for different species of 
inner crust matter in neutron stars, taking basic features 
of SQM and normal matter and providing reaction rates 
between them. No assumptions are made about the rela- 
tive amount of strange quark matter present in the inner 
crust, beyond postulating its existence; the pycnonuclear 
reaction calculations are limited to a single primitive lat- 
tice cell, with a single strangelet surrounded by normal 
matter ions. Future work may leverage the magnitude of 
the differences in reaction rates to modify solar heating 
models for neutron stars, perhaps providing clarification 
on the existence of strange quark matter. 

The paper is organized as follows. Section |TT] pro- 
vides an introduction into nuclear reaction rates in high- 
density environments. The characteristics of strange 
quark matter that appear to affect the pycnonuclear re- 
action rates are presented in Section Hill This is followed 
by a review of observed trends and an exploration of pos- 
sible future work presented in Section ITVl 



II. PYCNONUCLEAR REACTION RATES 

Generally speaking, in pycnonuclear reaction rates, the 
incoming and target particle densities usually considered 
in a nuclear reaction are replaced by geometric crystal 
lattice calculations that account for the separation dis- 
tance between a specific number of nearby, bound ions. 
The kinetic energy, usually involving a calculation of the 
relative velocity, is subsumed in an analysis of the total 
potential energy of a reacting pair of particles, stemming 
from both the configuration of the lattice and the vibra- 
tion of the bound ions. In turn, the incident flux through 



a pure coulomb barrier is replaced by an approximate 
calculation of the wave function transmitted through the 
total, screened, effective barrier between the two ions. 

In a high-density environment, such as a neutron star, 
it is expected that matter will settle into a low-energy 
configuration, such as the body-centered cubic which is 
assumed to have the greatest binding energy per nucleus 
[3^. The repulsive framework of similarly-charged ions 
enveloping two reacting ions in a bcc lattice reduces the 
Coulomb barrier between them, resulting in a barrier po- 
tential that is suppressed in both magnitude and range 
when compared to that denoted by a pure Coulomb law 
(4^ . Since ions trapped in such a lattice will continue 
to fluctuate about their equilibrium positions in their 
ground state, the ions may be able to penetrate the 
Coulomb barrier of a nearby ion (39j . In this effort, con- 
sistent with the most common treatment of pycnonuclear 
reactions within the model of the bcc lattice, the "static 
lattice" approximation - in which all the surrounding ions 
as well as the center of mass of the reacting ion pair are 
considered frozen at their equilibrium positions - is used 
pol . I4H . Furthermore, the electrons are assumed to act 
as a background electron gas with a density consistent 
for charge neutrality, moving within the arrangement of 
positive charge. 

Addressing the units of the reaction rate calcula- 
tion, lengths and energies for one-component (OCP) and 
multi-component plasmas (MCP) are usually measured 
in terms of the characteristic quantities r* and E* , re- 
spectively. These terms are deflned as 
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in analogy with the Bohr radius and the ground state 
energy of an electron in the Bohr model for hydrogen 
(40| . Here, M is twice the reduced mass for a pair of 
nuclei, and R is the atomic mass unit in grams. The 
density of the plasma can then be expressed in terms 
of r* and a dimensionless inverse-length parameter A, 
introduced by Salpeter and Van Horn in their seminal 
work on pycnonuclear reactions and defined as 
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where (Z) and (A) are the mean charge and mean mass 
number [3 Ell. 

With respect to the density model that describes the 
nuclear densities, a two-parameter Fermi (2pF) distribu- 
tion, or Woods-Saxon potential, is used: 



p{r) = 



Po 



l + cxp(r^) 
The normalization condition is 
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p (r) r^dr — X, 
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where X could be the number of protons, neutrons, or 
nucleons used to determine po [45|. This distribution 
can recreate the saturation of the nuclear medium and 
the rapid fall-off that brings out the notion of the radius 
Rq of the nucleus [45 1. 

On the whole, as an input to the 2pF distributions, the 
radii of most nuclei are well described by 



radius (in fermi) = 1.31 x A'^^^ - 0.: 



(9) 



where A is the number of nucleons in a given nucleus [46| . 
The rapid decrease in density is tied to the diffuseness pa- 
rameter a. This parameter plays a key role in subsequent 
models that vary the range of the nuclear force. Both 
nucleon and matter densities, coupled with their average 
diffuseness parameters Onucioon = 0.50 fm and amattcr = 
0.56 fm, respectively, give similar radii [47j . 



A. Long-Range Interactions 

The constituent calculations of a reaction rate are bro- 
ken up by where specific forces are dominant, with this 
factoring into long-range and short-range interactions 
possible because the lattice separations are much greater 
than the radii of the ions. As mentioned earlier, calcu- 
lations for the total kinetic energy involved in normal 
thermonuclear reactions are replaced by a calculation of 
the energies of the bound ions in pycnonuclear reactions. 
Beginning with the electrostatic lattice energy that main- 
tains the structure of the ions, an exact expression for the 
electrostatic screening potential is available, involving a 
six-fold sum of the locations of all nearest-neighbor nu- 
clei in the lattice (4^ . However, an approximation based 
on the static lattice paradigm is also possible, reducing 
the six-fold sum to a triple lattice sum based on the sep- 
aration distance r between the reacting ion pair. This 



approximation leads to the total electrostatic interaction 
energy per nucleus being given by 



bcc 

Coulomb 



z\^ 

= 1.81962 = 1.81962 A^;*, (10) 



where a is a constant of the lattice derived from the 
baryon number density [44| . However, the lattice en- 
ergy is not the only energy to consider; the nuclei are 
also vibrating in their ground states within the bcc lat- 
tice. The oscillation frequency of the nuclei around their 
equilibrium lattice sites is of the order of magnitude of 
the ion-plasma frequency ujp, defined by [40l | 



2.7426 X^^^E, 



Coulomb 1 



(11) 



with an anisotropic harmonic oscillator potential describ- 
ing the relative motion of the two ions. The energy- 
level spacings for vibrations in the z-axis separating two 
nearest-neighbor ions in the bcc lattice are |40j 



huj^ — 0.6752 hujp. 



(12) 



This vibrational energy in the z-direction along the axis 
connecting the two reacting nuclei can then be added 
to the electrostatic energy per particle in the bcc lattice 
to determine the total energy of the pair-wise nuclear 
reaction. 

To determine the wavefunction of the incoming parti- 
cle using the total, screened Coulomb potential, Salpeter 
and Van Horn used a Wentzel-Kramers-Brillouin (WKB) 
approximation to show that 



^(r) = 0.553 
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X exp 



with 



J~ X^^^K = 2.638- 3.6 A^/^ 
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is the wave function evaluated at the nuclear radius for 
the relative motion of the reacting nuclei in the static 
approximation of the lattice At this point, all of the 
long-range interactions for the pair-wise reaction rate W 
are included, resulting in 
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From here, only the eight nearest-neighbors surround- 
ing any given lattice point are taken into account, result- 
ing in the following expression for the reaction rate per 
unit volume: 
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Possible interactions with all other nuclei are ignored. 
The temperature-independent pycnonuclear reaction rate 
for a bee lattice, encompassing all the long-range effects 
of a pycnonuclear reaction, is thus given by 
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exp 



A1+A2 
-2.638 
Ai/2 



f3.90 X 10"*^ s 



(17) 



where p is in g-cm ^ and the S-factor is in units of 
MeV-barns [10, Hi. 



B. Astrophysical S-Factor 

The calculations up to this point have started from a 
macroscopic visualization of the long-range effects. The 
calculations that follow will invert this viewpoint, in- 
stead constructing the astrophysical S-factor from sep- 
arate models of the different forces with short-range ef- 
fects. With respect to units in S-factor calculations, the 
r* length scale of the bcc lattice is dropped for the more 
convenient use of fermi. One begins with the total ef- 
fective potential, the sum of the nuclear and unscreened 
Coulomb components [4ll |: 



^eff (r, E) = Vk (r, E) + Vc (r) + 



2^ 



(18) 



Since the dominant interaction at the low energies of py- 
cnonuclear reactions is via spherically-symmetric s-wave 
scattering, I is set to zero in all calculations presented in 
this effort 0. 

The first term, the nuclear potential, is based on one of 
several methods created in a recent model, referred to in 
this work as the Sao Paulo model, which was developed 
by a joint collaboration from several universities and in- 
stitutes in Brazil [45, 48, 49,, ,5fl,]. One starts with the 
concept that the effective one-body interaction between 
two nuclei can be written schematically as 

l/(r,r") = ybare(r,r-') 

+ Y,V,{r)G,{r,r':E)V,{r'). (19) 

i 

The first term on the right is the bare interaction, rep- 
resenting the ground state of an interaction between two 
nuclei. The second highly-energetic term encompasses 
polarization and intermediate states, like inelastic chan- 
nels [4^. This second term is not considered in this effort. 

The first term is the focal point of the Sao Paulo model 
and can be calculated using multiple methods, though 
the one presented here associates an energy dependence 
of the derived nuclear potential with nonlocal quantum 
effects. It begins by concentrating on the effects arising 
from the Fermi nature of the nucleons. When calculat- 
ing interaction potentials between nuclei, these effects 



translate into a nonlocality. The nonlocality of the bare 
interaction is solely due to the Pauli exclusion principle, 
hereafter referred to as Pauli nonlocality [i^. This Sao 
Paulo methodology creates a local equivalent potential 
that makes use of a double-folding methodology whereby 
a representation of the effective nucleon-nucleon force is 
folded, i.e., integrated, with the densities of the interact- 
ing nuclei: 



old 1: 



Pi {ri) 92 {t2) "nn drx dr2 ■ (20) 



The result is dependent only on the number of nucleons 
in the nuclei. 

The Sao Paulo group adopts the following form for the 
bare nuclear potential between two nuclei: 



Vbare (r, r') = V^NL ( 7, ' ] ^Xp 



{r-r'Y 

62 



(21) 



where h is the range of the Pauli nonlocality of the ion-ion 
interaction 48]. In a sense, this term measures the range 
of the nonlocal effect where the nuclear force is operative. 
Jackson and Johnson explain in previous work that 



60 m 



(22) 



within their single folding model, where 60 is the nucleon- 
nucleus nonlocality parameter, m is the nucleon mass, 
and p is the reduced mass of the nucleus- nucleus system 
[sH - Fitting to nucleon-nucleus scattering data, Percy 
and Buck later show that 60 = 0.85 fm '521. Contin- 
uing to use the density-folding formalism, the nonlocal 
potential Vnl in Eq- (ED) is written as [1^ 



Vnl = Vfoiding = J {fi)P2 {r2) VNN dri dr2 ■ (23) 

Using the optical model, the bare nuclear potential is 
rewritten in terms of a central potential. This requires 
extracting a local-equivalent potential from a nonlocal 
expression. A local equivalent potential is defined as (s^ l 



V (f, r') V^E {r') dr' = V (f, E) V-e {r^ ■ 



(24) 



Accordingly, in the central potential notation of the opti- 
cal model, the total nuclear interaction can be rewritten 
as 'A7i 



1 r°° 

— Vi{R, R') ui {R') = Vle (R, E) 
\R) Jo 

+ iWLEiR,E) . (25) 



ui 



As can be seen, there is an energy-dependent, local equiv- 
alent potential in the above equation for every term in the 
effective one-body interaction, with the first term being 
predominantly real and equivalent to the nonlocal, bare, 
nuclear interaction while the second term is complex (and 
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still ignored). In this model, the energy dependence of 
Vle {R, E) is chiefly due to Paiili nonlocality in the bare 
nuclear interaction Vbaic {r, r'). 

To create this local equivalent potential, the Sao Paulo 
group uses the Percy and Buck prescription, namely that 
for a given non-local optical potential depth (V^nl), the 
corresponding local potential depth (Vle) at each energy 
should be chosen to satisfy the approximate relation 



Vle exp 



(26) 



where M is the mass [52] ■ Using Eq. (I21|) and recognizing 
that T4viL determines the nonlocal potential depth for the 
bare nuclear interaction, the local equivalent expression 
in this case can be written as the following [isj : 



VLEiR,E) « FF(i?)e-T[^-^'=(^)-^^^(^^^)l, (27) 

2n^' 
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(28) 



This expression has multiple possible interpretations. 
In one interpretation, the exponent can be tied classically 
to the kinetic energy E and relative speed v between the 
nuclei by 



resulting in 



'-Ek (R) ^-[E-Vc (i?) - Vle {R, E)] , (29) 



Vle {R, E) 



V-p (R) e-[™«''oi'/2?i]= 



(30) 



where c is the speed of light. This interpretation equiv- 
alently relates the effect of Pauli nonlocality to that of 
a velocity-dependent nuclear interaction. In the current 
work, the local-equivalent potential is instead related di- 
rectly to the folding potential, which incorporates an ef- 
fective nucleon-nucleon interaction dependent on the rel- 
ative speed V between the nucleons [47] , 



Vle {R,E)^Vy 



Vf = 



Pi {ri) P2 (^2) WNN [v, R-fl+ fij dfl (1^2 



WNN (w, r)^Vi (f) e ""'/"^ 



(31) 

(32) 
(33) 



This additional modularity allows for different nucleon- 
nucleon interaction models to be incorporated in future 
work. 

Accordingly, Vf (f) is the next term to be considered. 
The Sao Paulo group extrapolated that the effective 
nucleon-nucleon interaction could be derived by reusing a 
folding potential from previous heavy-ion potential anal- 
yses, 



Vf ■ 



Pm (ri) Vq d(^v,R~fl+ r^^ dfl df^, (34) 



where Vq = —456 MeV fm . In this scheme, pm is the 
matter density of the nucleon, for which they assumed 
an exponential shape. 



Pm (r) = po e-'-/"" 



(35) 



based on electron scattering experiments that determined 
the intrinsic charge distribution of the proton in free 
space. Again, po can be determined by the normaliza- 
tion condition from the density calculations 45] . The in- 
tegration of Vf results in a finite-range, nucleon-nucleon 
interaction. 



Vf ■ 



Vo 
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— + ) , (36) 



shown in Fig. (TJ where the best-fit value for the matter 
diffuseness Om of a nucleon is 0.3 fm. 

When multiplied by the velocity-dependent exponen- 
tial in Eq. p3|) . effects due to Pauli nonlocality are in- 
corporated into this local equivalent, nucleon-nucleon in- 
teraction. Unfortunately, with its focus on recreating 
nonlocal effects, this methodology does not recreate a 
repulsive core for the nucleon potential, an empirical fea- 
ture that gives rise to the saturation of the nuclear force. 
However, this method provides a good fit of heavy-ion 
potential strengths and is similar to the M3Y potential 
in the surface region (isl ]. 
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FIG. 1: Nucleon potential without Pauli nonlocality as a func- 
tion of distance from nuclear center. 

The nuclear potential of a system of two C12 ions can 
be seen in Fig. [2l The target ion is considered locked in 
place near the far-right edge of the graph. The distance 
between target and projectile centers at the right edge of 
the graph is greater than the actual distance between the 
lattice points, i.e., in a two-body picture, the projectile 
has actually moved past the target; this is the cause of 
the slight uptick in the nuclear potential. 

Moving on from the complications of the nuclear po- 
tential, the short-range effects of the Coulomb interaction 
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FIG. 2: Nuclear potential between two C12 ions as a func- 
tion of position on the r-axis between nuclear centers for 
p = 10^^ g cm-^ 



in the effective potential are taken into account. Outside 
of the nuclei, the point-particle assumption holds. How- 
ever, in the interior of the nucleus, the expression must 
be modified to account for the finite size of the nucleus, 
and the potential energy due to the Coulomb repulsion 
becomes [43 
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(37) 



within the nucleus. 

Having incorporated nuclear and Coulomb interactions 
into a total effective potential, the calculation of the S- 
factor is accomplished through the use of a barrier pen- 
etration model. This heuristic treatment begins with a 
calculation of the transmission coefficient through the to- 
tal effective potential, incorporating this into a model 
based on fits of thermonuclear data to compute the nu- 
clear cross section j41|]. The definition of the S-factor, 



S{E) = a{E)xEx 



(38) 



with the Gamow coefficient T defined for transmission 
through an unscreened Coulomb barrier, can then be 
used to calculate the reaction rate. 

The transmission coefficient is derived by way of a 
WKB approximation that requires as input the total ef- 
fective potential and the total energy. Using the momen- 
tum of the reacting system, possible trajectories through 
the total potential are orthogonal to the surfaces of con- 
stant phase of the wavefunction [ssj. To solve for the 
phase, the momentum is integrated between the classical 
turning points of the effective potential, using 



k{x) - ^{2m [Vesix) - Eo]}'^^ 



(39) 



as the momentum in one dimension [53 
dimensional calculation for phase is then 



The three- 



WKBi = ±jj %s{x)-Eo] dr. (40) 

From this, the probability corresponding to the passage 
of the particle through the effective potential is given by 
the transmission coefficient 



T, = 



IV'tr 



1 + e 



WKBi 



(41) 



In accordance with the barrier penetration model, the 
fusion cross section is tied to the summation of the to- 
tal particle flux transmitted through the barrier up to 
the greatest value of angular momentum in the effective 
potential [4ll |. 



a{E) 
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l)Ti. 



(42) 



The astrophysical S-factor is then calculated by multiply- 
ing a (E) , the total energy, and the inverse of the trans- 
mission coefficient through a pure, unscreened Coulomb 
barrier. 

The fact that the inverse of the Gamow penetration 
factor is being used is worthy of a second look. The cal- 
culations for the total cross section of the reaction have 
explicitly included a pure Coulomb interaction in the to- 
tal effective potential. It may seem as if these efforts 
are being undone by removing the Gamow penetration 
factor from the S-factor, but this is not the case. The 
inclusion of the Coulomb interaction in the effective po- 
tential was necessary to find the correct turning points 
for the WKB approximation needed to calculate the to- 
tal cross section, particularly the turning point for the 
onset of the nuclear force. But, it is important to re- 
member that the S-factor is supposed to be an expres- 
sion that compartmentalizes nuclear force effects. Di- 
vesting the S-factor of Coulomb forces outside of the nu- 
clear radius follows this paradigm. Furthermore, the re- 
action rate calculation separately accounts for the effects 
of the more accurate, screened Coulomb potential of the 
bcc lattice. It is necessary to divide out the traditional, 
unscreened Gamow penetration factor in order to avoid 
double-counting Coulomb effects. 



C. Atomic Matter Pycnonuclear Reaction Rates 

With all of the constituent parts for the calculation in 
place, pycnonuclear reaction rates for various elements 
are provided in Figs. [2]and[3]for a range of densities. The 
reaction rate results for OCP and MCP nuclear matter 
are a good match to those calculated by the originators 
of the model, an extended collaboration that included 
the Sao Paulo group that created the nuclear interaction 



7 



methodology and the Joint Institute of Nuclear Astro- 
physics ( JINA) at the University of Notre Dame [43|] . 

At this interim point, the pycnonuclear reaction rates 
confirm an expected trend. With the data, one can begin 
to see that the reaction rates decrease as the atomic num- 
bers of the nuclei involved increase. When this method- 
ology is extended to heavy ions that may be found in the 
inner crust of a neutron star - potentially iron all the 
way up to neutron drip - the reaction rates continue to 
plummet, as can be seen by the rapid change in scale for 
Fig. 21 This makes sense: The increase in the Coulomb 
repulsion between larger nuclei would make it more difh- 
cult to tunnel through the barrier between the ions, even 
with the effects of Coulomb screening. 



40 



-80, 






+ 






— - ^^c 


+ 








+ 


^^0 





—I I l_ 



9.5 



10 



—I I I I I I i_ 



10.5 



11 



log^g Density (g/cm ) 



FIG. 3: (Color online) Reaction rates for OCP and MCP 
atomic matter as a function of density. 
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FIG. 4: (Color online) Reaction rates for OCP neutron-rich 
nuclei as a function of density. 



III. STRANGE QUARK MATTER 

From here, the effort expands to an attempt at the cal- 
culation of pycnonuclear reaction rates between normal 
nuclear matter and a simple approximation of strange 
quark matter. It is important to acknowledge immedi- 
ately that this extension to strange quark matter makes 
no attempt to modify the methodology of the strong 
interaction as it has been presented for atomic matter. 
However, several other models implemented for normal 
nuclear matter are modified to accommodate theoreti- 
cal characteristics of strange quark matter, particularly 
the bulk charge of a strangelet as well as its radius, den- 
sity, and baryon number. For the purposes of this effort, 
the term strangelet is used to refer to an approximately 
nucleus-sized strange quark matter nugget that replaces 
one of the two normal matter nuclei in the previous re- 
action calculations. 



A. Strangelet Mass 

To replace an atomic matter nucleus within the bcc 
lattice with a strangelet, one needs to know the approx- 
imate mass of the strange quark matter nugget. There 
are numerical studies that have modeled the late stages 
of in-spiral in binary systems of neutron stars and mas- 
sive objects such as black holes, white dwarf stars, or 
other neutron stars. As tidal forces and, eventually, the 
actual collision overcome the surface tension of strange 
quark matter and rip the compact objects apart, the re- 
lease of macroscopic lumps of matter between 10~* and 
10""'^ solar masses, or A oc 10'^^, seems to be a charac- 
teristic feature [54|- It seems a reasonable assumption 
that a significant fraction of this tidally-released mate- 
rial remains trapped in orbit with typical speeds around 
0.1c and subject to frequent collisions. If one calculates 
that the kinetic energy from these collisions is converted 
to supplying the extra surface and curvature energy nec- 
essary for forming smaller strangelets, then the released 
mass may end up in strangelets with A w 10^ — 10^ [s^ . 
This range is the focus of the input parameters for strange 
quark matter baryon numbers in the pycnonuclear reac- 
tion rate calculations presented below. Since strangelets 
are self-bound, their mass-radius relationship is given by 
M = edV 47rei?3/3 Also, instead of using the 
more rigorous two-parameter Fermi model, the density 
is assumed to be approximately constant throughout the 
strangelet. 



B. Strangelet Charge 

The overall charge of the strangelet is the final parame- 
ter changed, and the one that most affects the outcome of 
the pycnonuclear reaction rates in this study - again, ac- 
knowledging that the nuclear interaction is unchanged. 
We assume that the strangelets are made up of either 
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ordinary strange quark matter or color superconducting 
strange quark matter whose condensation pattern is the 
color-flavor-locked (CFL) phase. Once crucial difference 
between non-CFL (NCFL) and CFL quark matter is the 
equality of all quark Fermi momenta in CFL quark mat- 
ter which leads to charge neutrality in bulk without any 
need for electrons (ssj . This has most important conse- 
quences for the charge-to-mass ratios of strangelets. For 
ordinary strangelets, the charge is approximately 



Z K 0.1 ml^o A, ^<10^ 



Zk8 



'150 



A > 10^ 



(43) 



(44) 



where rriiso = nis/ 150 MeV and rris is the mass of the 
strange quark. For small A, the charge is the volume 
quark charge density multiplied by the strangelet vol- 
ume with a result that is proportional to A itself. This 
relation holds until the system grows larger than around 
5 fm, or A ^ 150, at which point the charge is mainly 
distributed near the strangelet surface, and Z cx A^/'^ 
[5^. In contrast to this, the charge-to-mass ratio of CFL 
strangelets is described by [s^ 



0.3 mi 50 A^/^ 



(45) 



Recently it has been shown that CFL strangelets may 
even have zero electric charge if the quark pairing is 
strong [53|. If strangelets should exist in the crusts of 
neutron stars, this difference may provide a test of color 
superconductivity. 



C. Reaction Rates with Strange Quark Matter 

After incorporating these changes due to strange quark 
matter, pycnonuclear reaction rates between normal nu- 
clear matter and matter nuggets that incorporate some 
characteristics of strange quark matter are compared in 
Figs. [5] to [H In Figs. [5] and [3 the atomic nuclei are 
paired with strangelets of equivalent baryon number; in 
Figs. Elandini strangelets of baryon number 500 are used. 
Previously-shown data for the reaction rates between two 
atomic matter nuclei are also provided in the graphs for 
comparison. In these calculations, the strange quark was 
assumed to have a mass of 300 MeV. 

Again looking for possible trends, one can see that the 
pycnonuclear reaction rates still decrease as the atomic 
numbers of the nuclei or strangelets involved increase. 
Additionally, the reaction rates between atomic nuclei 
and the strange quark matter nuggets are very different 
than that of normal nuclear matter alone. However, the 
possibility of a color superconducting state for strange 
quark matter has dramatic consequences. The extent of 
the increase in reaction rates when compared to those 
of atomic nuclei alone is noteworthy, with escalations of 
over twenty orders of magnitude, at a minimum, at most 
densities of interest in the cases studied. By comparison. 
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FIG. 5: (Color online) Reaction rates for ^ Fe and a strange 
quark matter nugget with 56 baryons as a function of density. 
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FIG. 6: (Color online) Reaction rates for ^®Fe and a strange 
quark matter nugget with 500 baryons as a function of density. 



non-CFL (NCFL) strange quark matter reaction rates 
are always significantly lower than those of CFL strange 
quark matter with equivalent baryon number, because of 
the electric charge difference. The combination of these 
divergent effects - with CFL and non-CFL strange quark 
matter rates maintaining a significant gap between them, 
the inclusion of CFL strange quark matter increasing the 
reaction rates of atomic nuclei, and increasing baryon 
number decreasing reaction rates for all systems - makes 
it difficult to determine a priori how reaction rates for a 
given system will compare to each other. 
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FIG. 7: (Color online) All reaction rates calculated for ^''Fe 
as a function of density. 
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FIG. 8: (Color online) Reaction rates for ^^*Kr and a strange 
quark matter nugget with 118 baryons as a function of density. 



D. Quark Nugget Concentrations in the Crusts of 
Neutron Stars 



An upper limit on the quark nugget concentration in 
the crustal regions of neutron stars accreting matter from 
white dwarfs would be set by the quark nugget flux lev- 
els that could have been left over from the Big Bang or 
from collisions of strange stars [2^, [2^, IH, [5^ . The cap- 
ture of these quark nuggets by main-sequence stars would 
be an inevitable consequence of the strange quark mat- 
ter hypothesis. Due to their large radii, main-sequence 
stars arise as large-surface long-integration-time detec- 
tors for the quark nugget flux EH- The accreted 
nuggets should be thermally distributed throughout the 
interior of main-sequence stars, leading to quark nugget 
concentrations that would monotonically grow over the 
lifetime of main-sequence stars. For a constant bath of 
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FIG. 9: (Color online) Reaction rates for ^^*Kr and a strange 
quark matter nugget with 500 baryons as a function of density. 



isotropic monoenergetic quark nuggets, the baryon num- 
ber accreted onto a main-sequence star can be estimated 
from [13, m 



A < 1.6 X 10' 



47 



M 

Mo 



1 + 0.164 ( — 



-0.15 



M 



"250 P24 
-0.25\ 



(46) 



where f25o = ^00/250 km/s and p24 = Poo/lO"^'' g/cm'^. 
The quantities Voc and denote the quark nugget speed 
and contribution to the density of the galactic halo far 
from the star, respectively, and M denotes the mass of a 
main-sequence star in units of the mass of the sun, Mq. 
If one assumes that all the dark matter in the halo of our 
galaxy would consist of quark nuggets, which is certainly 
an overestimate, then P24i'25o ~ 1- Equation then 
leads for typical progenitor star masses to A ^ 10'*^ for 
the total baryon number of the quark nuggets in white 
dwarfs, formed from main-sequence stars. This value is 
to be compared with the typical baryon number of the 
crust of a neutron star, which is around 10^^. Hence, 
based on this scenario, one would expect less than one 
quark nugget per 10^ atomic nuclei. Of course, this ratio 
may be significantly greater in regions of the Universe 
(e.g. Globular Clusters) where collisions of strange stars 
would have generated quark-nugget flux levels that are 
greatly enhanced in comparison to the homogeneously 
spread quark-nugget flux level assumed above. 



IV. CONCLUSIONS 

The methodology presented in this paper provides in- 
sight into the magnitude of the rates of pycnonuclear 
reactions between a simple model of strange quark mat- 
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ter and atomic nuclei. Recreating a nuclear interaction 
model and incorporating it in high-density nuclear re- 
action processes, this effort reconstructs previous calcu- 
lations for normal nuclear matter. It then expands in 
scope, extending the model to estimate the changes in 
the reaction rates as distinctive theoretical characteris- 
tics of strange quark matter are substituted for atomic 
matter. 

Assessing the results for all of the data presented, there 
is an unsurprising decrease in the reaction rates for nuclei 
with larger atomic numbers, corresponding to the ampli- 
fication of the Coulomb barrier as the positive charge 
is increased. In the extended model, though the charge 
suppression of strangelets was expected to have an effect 
on pycnonuclear reaction rates, the increase by at least 
twenty orders of magnitude for reactions between CFL 
strange quark matter and atomic nuclei at the densities 
of a neutron star crust was surprising. Reactions that 
incorporated non-CFL strangelets maintained rates that 
were significantly lower than the CFL results, another in- 
dicator of the additional reduction in charge for the color 
superconducting system. With these conflicting trends 
and the associated inability to presciently determine the 
reaction rates of the differing systems with respect to 
each other, this methodology for calculating pycnonu- 
clear reaction rates will be of great importance for future 
astrophysical studies. 

As an example, we mention the incorporation of the 
fusion reaction rates computed in this paper into pyc- 
nonuclear reaction fusion network calculations, such as 
performed by JIN A [6^. Such calculations will pro- 
vide key information about the heat release from neu- 
tron stars whose crusts may be contaminated with (su- 
perconducting) strange quark matter nuggets. Since the 
reaction rates of the latter differ by many orders of mag- 
nitude from the reaction rates of ordinary nuclear crust 
matter, nuclear fusion network calculations may have 
the potential to shed light on the existence of strange 
quark matter nuggets and on whether or not they are 
in a color superconducting state, as suggested by QCD 



A. Uncertainties 

There are several areas of uncertainty. First, the py- 
cnonuclear reaction rate equation in Section |TT] makes 
use of the WKB integrals J and K, integrals that 
were numerically calculated for densities between 1.7 x 
IQ-^g cm-3 < p/{A)A^Z^ < 1.7 x lO^g cm-^ The 
densities possible from the outer crust to the core of a 
neutron star are all several orders of magnitude greater 
than this density range. In the attempt to recreate previ- 
ous efforts that used and did not change the pycnonuclear 
reaction rate equation, the applicability of the Salpeter 
and Van Horn calculation at these densities has been as- 
sumed. 

There is also the issue of the nuclear potential model 



for atomic matter. Referring again to Fig. [Tl it should 
be noted that the Sao Paulo model does not imple- 
ment a repulsive, hard core for the nuclear interaction. 
Given its focus on recreating nonlocal effects with a local- 
equivalent potential, it is quite possible that it was never 
intended to do so. Though it shows similarities to the 
M3Y potential in the surface region of nuclei and pro- 
vides a good fit to data of heavy-ion potential strengths 
[4^ , its validity has been called into question |63, IBll . 
It is also critical to remember that this study leaves the 
nuclear interaction model unchanged when the adjust- 
ments for strange quark matter are made. Any possible 
conclusions must be weighed against this deficiency. 



B. Future Work 

There are several avenues available for future work. 
The initial suggestions begin with the known sources of 
error just mentioned in Sect. IIV Al Though any analy- 
sis and modification of the / and K WKB integrals for 
the pycnonuclear reaction rate equation is a task heavily 
steeped in numerical and complex analysis, Salpeter and 
Van Horn do offer many details about their calculations 
in their published work [13, [g^l ■ A verification that the 
current reaction rate equation properly accounts for the 
densities under consideration in a neutron star would be 
of use. 

There is another opportunity in that the modularity 
of this methodology allows for comparisons of different 
nuclear force models. If one does not intend to model ex- 
otic matter, the current methodology may still be used 
to assess changes in reaction rates due to various nuclear 
interaction models. It may be possible that such com- 
parisons may provide support for the use of one nuclear 
model over another. 

Future efforts may include an improved model of the 
nuclear interaction with strange quark matter. The cur- 
rent double-folding process that determines the total nu- 
clear potential uses a model of the nuclear force that is 
assumed to be applicable to both the target and projec- 
tile nuclei within the lattice. If a strangelet interaction 
model is included, future efforts may need to alter the 
subordinate nuclear potential vnn prior to its integra- 
tion with the strangelet and atomic matter densities. On 
the other hand, a simpler analysis could focus on any ef- 
fects of the reaction rates that may stem from changing 
the mass of a strange quark, assumed to be 300 MeV, 
within the current strange quark matter model. 

Separately, it may also be useful to assess the effect of 
strangelets or different atomic nuclei on the pycnonuclear 
reaction rates of nearby interacting particle pairs. This 
sort of calculation, focused on the effects of lattice imper- 
fections, has not been extensively studied even for nor- 
mal nuclear matter. This is not a three-body calculation, 
since the environment already contains many nearby par- 
ticles frozen in a lattice. The lattice imperfection would 
first alter the geometry of the electrostatic screening po- 
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tential in Section |TT1 affecting the vibrational energy lev- 
els of a nearby interacting ion pair and, in turn, influenc- 
ing the reaction rate. Since the screened Coulomb po- 
tential is also used to determine the wavefunction of the 
incoming particle, it is likely that this particular line of 
research would require extensive changes to the pycnonu- 
clear reaction rate equation itself. The methodology for 
the calculation of the S-factor could remain unchanged. 

In addition, follow-on efforts may wish to incorporate 
pycnonuclear reaction rates into stellar thermal models, 
possibly flowing potential characteristics of strange quark 
matter to observable predictions. As a first approxi- 
mation, the zero-temperature pycnonuclear reaction rate 
was used in this effort; future efforts may wish to incor- 
porate the temperature-dependent pycnonuclear reaction 
rate equation developed by Salpeter and Van Horn [ioj . 
One could build on this by calculating an energy genera- 
tion rate in units of erg g~^ s"'^ that can be obtained by 
multiplying the pycnonuclear reaction rate equation by 
Q/2p [40]. Q is the energy release in a single reaction. 
Approximations for Q would need to be made for strange 
quark matter, perhaps building on estimates that the to- 



tal energy release for normal nuclear matter is 1.45 MeV 
per baryon [6^1 . These thermal calculations could then 
be incorporated into stellar heating models, evaluating 
possible variations in total heat output for neutron stars 
due to normal matter or strangelets that may be present 
as a result of state changes or accretion. Differences in 
reaction rates may thus provide a clue to the presence of 
strange quark matter in astrophysical environments. 

In any case, it is hoped that this investigation may 
serve as a stepping stone for further studies, and that 
these future possibilities may tie astrophysical obser- 
vations of high-density objects to speculations about 
the possible existence of (color superconducting) strange 
quark matter. 
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